Subgap transport in ferromagnet-superconductor junctions due to magnon-assisted 

Andreev reflection 
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We propose a new process of magnon-assisted Andreev reflection at a ferromagnetic metal 
- superconductor interface, which consists of the simultaneous injection of a Cooper pair from 
the superconductor and the emission of a magnon inside the ferromagnet. At low temperature 
this process represents an additional channel for sub-gap transport across an FS interface, which 
lifts restrictions on the current resulting from the necessity to match spin-polarized current in the 
ferromagnet with spin-less current in the superconductor. For a junction between a superconductor 
and a ferromagnet with an arbitrary degree of polarization, the inelastic magnon-assisted Andreev 
reflection process would manifest itself as a nonlinear addition to the I{V) characteristics which is 
asymmetric with respect to the sign of the bias voltage and is related to the density of states of 
magnons in the ferromagnet. Expressions for the subgap I{V) characteristics are given for arbitrary 
interfacial quality whilst the limiting cases of uniformly transparent and disordered interfaces are 
discussed in detail. 

PACS numbers: 72.25.-b, 75.30.Ds, 74.80.Dm. 



I. INTRODUCTION 

Spin polarized transport is a subject of intense re- 
search, motivated by the desire to develop a form of elep.- 
tronics which utilizes the spin polarization of carrierso. 
Ferromagnetic (F) metals have more carriers of one spin 
polarization (known as majority carriers) present at the 
Fermi energy Ep than of the inverse polarization (mi- 
nority carriers) and, in a ferromagnet-ferromagnet (FF) 
junction, the resistance depends on the relative orien- 
tation of the magnetization in the ferromagnets. With 
parallel magnetizations, carriers flow from majority to 
majority bands (and minority to minority) whereas in a 
junction with antiparallel magnetizations, carriers flow 
from majority to minority bands (and vice versa). The 
resulting spin current mismatch produces a larger contact 
resistance in the antiparallel orientation, an effeclj-known 
as tunneling magnetoresistance in FF junctionals and 
giant magnetoresistance in multilayer structuresoE . 

Meanwhile, recent achievements in fabrication tech- 
niques have led to the possibility of creating hybrid het- 
erostructures combining ferromagnetic (F) and supe|rcon- 
ducting (S) elements and it has been pointed outtlia that 
the effect of spin current mismatch may affect thefCoa- 
ductance of a ferromagnet-superconductor junctionLlLj. 
At low temperatures and small bias voltage (subgap 
regime, T ^ A, eV < A), curreut, flows through the 
interface due to Andreev reflectionlla whereby particles 
in the ferromagnetic region with excitation energies e 
smaller than the superconducting gap energy A are re- 
flected from the interface as holes. Since subgap trans- 
port in the superconductor (S) is mediated by spinless 
Cooper pairs the spin current is zero in the supercon- 
ductor in contrast to the ferromagnet. In high-quality 
metallic type FS junctions, this has been shown to result 



jn accumulation in the vicinity 
Such an accumulation of non- 



in a non-equilibrii 
of the FS interfaci 
equilibrium spin density generates a compensating spin- 
current flow and, therefore, an additional non-local in- 
terfacial resistanceQ foaaiied at the spin-relaxation length 
scale in a ferromagnemO. 

An increasing number of experiments have been de- 
voted to studying the transport properties of feccomag- 
net/insulator/superconductor junctions. SomeOu have 
measured the degree |0£j5pin-polarization of various fer- 
romagnets and otherscltJ have studied the effect of spin- 
polarized quasiparticle injection into high- Tc supercon- 
ductors. In tunnel FS junctions, the necessity to match 
microscopic spin currents at the interface creates an op- 
portunity for complex electron transfer processes to man- 
ifest themselves in the currenl formation. For instance, 
spin relaxation phenomenaoo, such as spin-orbit scat- 
tering at impurities or magnon emission, can reduce the 
spin current mismatch. In a given junction, spin-orbit 
scattering (which is an elastic process) would cause a 
linear recLurtinn in the value of the additional contact 
resistanceED'EJ. On the other hand the inelastic process 
of magnon emission would manifest itself as a modifica- 
tion of the form of the I{V) characteristics. Indeed, non- 
linear I(V) characteristics due to magnon-assisted tun- 
neling between two feipoHiagnets have already befp-stud- 
ied both theoreticallj£3~EJ and experimentallyEj'E£l with 
a view to relate the second derivative of the current to the 
density of states of magnons n{uj) as (Pl/dV^ oc n{eV). 

The aim of this paper is to investigate Andreev reflec- 
tion accompanied by the emission of magnons in tunnel 
contacts between a ferromagnet and a superconductor (a 
preliminary discussion of the half-metallic ferromagnet- 
superconductor case is given in Ref . . For this purpose 
we calculate the nonlinear subgap current / using the 
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tunneling Hamiltonian method i3 and the nonequilibrium 
Green functions techniqud^Ju generalized for describing 
contacts with ferromagnetic electrodes. The differential 
conductance G{V) — dl/dV {V is the voltage) is shown 
to be the sum of the usual Andreev conductance Ga and 
a contribution Gin from inelastic processes in the ferro- 
magnet: G = Ga + Gm- For a uniform partially tr|a 
parent interface Ga is given at low temperatures bj; 



W « 1, 



(1) 



where the factor stands for the probability of the 
simultaneous tunneling of two particles involved in the 
process of Andreev reflection, h is the Planck constant 
and e is the absolute value of the elementary charge. 
Throughout this paper we use the notation f and J, to 
represent majority and minority spin bands, respectively, 
and the factor in Eq. (|l|) denotes the area of the max- 
imal cross section of the Fermi surface of minority (spin 
'down') electrons in the plane parallel to the interface, so 
that Andreev conductance is proportional to the number 
of minority spin channels AIli jh? \a. a. contact with area 
A. 

The inelastic contribution Gin to the differential con- 
ductance results from the interaction between spins of 
the conduction electrons and the quasilocalized moments 
of the inner atomic shells which are responsible for the 
magnetism. In the second quantization language emerg- 
ing from the Holstein-Primakoff transformation 3Lthis in- 
teraction - treated within the s-f (s-d) modeO'EEl ~ rep- 
resents magnon emission or absorption accompanied by 
a spin-flip process which is preceded or followed by An- 
dreev reflection. In the zero-temperature limit where the 
magnon emission is only allowed the inelastic conduc- 
tance Gin is given for a uniform interface by 



G^n{V) = —\t\' 



M{eV)^ + Mi-eV)^ 



2x 



M{x) = /\^Q{x) I n{Lu)duj X 
Jo 

1 

X 



(2) 
(3) 



A2-a;2 A^-{x-i 



2eV < 2A < Es- 

The first term in Eq. represents the contribution from 
the majority (spin 'up') states close to the Fermi level 
and Hi denotes the maximal cross section of the major- 
ity Fermi surface in the plane parallel to the interface. 
The second contribution is related to electronic states lo- 
cated in the vicinity of the minority Fermi surface. It 
is assumed that the magnon local density of states near 
the interface, 0(w) in Eq. (^, is non-zero for tu > luq, 
where luq is a relativistic anisotropy gap in the magnon 
spectrum which may arise from a spin-orbit coupling or 
from the demagnetization geometry of an F-layer (see. 



e.g., Ref ^m. S is the spin of the localized moments. 
In Sectioii|| we describe in more detail the model and 
technique used, and we give general results for the I{V) 
characteristics valid for arbitrary interfacial quality. 

Both terms in the square brackets in Eq. (|^) are asym- 
metric functions of voltage (this is emphasized by the 
presence of the unit step function 8 in Eq. (||)). The 
first term is non-zero for positive voltage V > LUo/(2e), 
while the second one is non-zero for negative voltage 
V < —uJo/(2e). This feature can be explained using the 
sketch in Fig. 1 which illustrates the tunneling process 
between an S electrode on the left hand side and an F 
electrode on the right for (a) ^ > and (b) V < 0. 
We have adopted the convention that positive (negative) 
voltage results in a Fermi energy Ep in the ferromagnet 
that is lower (higher) by energy \eV\ than the Fermi en- 
ergy Es in the superconductor. For V > 0, Fig. 1(a), 
Andreev reflection results in the injection of both a ma- 
jority (spin 'up') and a minority (spin 'down') electron 
into the ferromagnet. The electron energies, ±e with re- 
spect to Es and eV ± e with respect to Ep, are both 
above Ep because of the need to move into unoccupied 
states in the ferromagnet, e < eV. Due to the existence 
of the intra-atomic exchange interaction with the local- 
ized moments, a dynamic process, shown schematically 
in Fig. 1(a), allows a spin down electron to tunnel into a 
virtual, intermediate spin down state above Ep and then 
emit a magnon, depicted in Fig. 1(a) as a fiip of a core 
spin, enabling the electron to incorporate itself into an 
empty state in the majority, spin up conduction band. 
Since the magnon carries away spin equal to -1, the total 
spin is equal to zero in the above process. The relevant 
magnon energies u are less than the energy of a Cooper 
pair 2eV measured with respect to the chemical poten- 
tial in the ferromagnet (see Eq. (^)). Since the super- 
conducting gap energy A is usually much less than the 
magnon Debye energy, the inequalities cu < 2eV < 2A 
imply that mainly long-wavelength magnons with mo- 
menta much smaller than the Fermi momentum assist 
the subgap transport. This means in turn that only elec- 
tron states near the majority Fermi surface are involved 
in the magnon-assisted transport at V > 0. 

On the other hand, for 1/ < 0, Fig. 1(b), Andreev re- 
flection leads to the injection of a spin up and a spin 
down electron from the ferromagnet into the supercon- 
ductor with energies e above and below Es- Due to 
the intra-atomic exchange the following dynamic process 
can contribute to current formation at zero temperature 
T = 0: one of the spin down electrons below Ep emits 
a long-wavelength magnon and forms a virtual, interme- 
diate spin up state. Such a process allows it to tunnel 
into the superconductor along with another spin down 
carrier. Their energies, — |ey| ± e with respect to Ep, 
are both below Ep because of the need to have initially 
occupied states in the ferromagnet. Since emission of 
a long-wavelength magnon does not change significantly 
the energy and momentum of the spin down electron, the 
initial and intermediate electron states are close to the 
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minority Fermi surface. We will point out later that both 
the majority and minority magnon-assisted currents are 
quadratic in the interaction constant F relevant to the 
intra-atomic exchange. On the other hand, they are also 
quadratic with respect to the electron lifetime in the vir- 
tual state which is of the order of fi/V. For this reason 
Eq. (||) for the magnon-assisted conductance does not 
contain the interaction constant. 



— [G{V)-G{-V)] 



Ae^A\t 




n(\2eV\) 



l-(|ey|/A) 
eV 



de 
2A 



(4) 
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FIG. 1 - Schematic representation of the magnon- 
assisted tunneling process between a superconductor on 
the left hand side and a conventional ferromagnet on the 
right for (a) 1/ > and (b) F < 0. For (a) F > a 
spin down electron tunneling into the ferromagnet may 
emit a magnon and incorporate itself into the majority 
(spin up) conduction band. For (b) < a spin down 
electron with energy below Ep may emit a magnon and 
form an intermediate spin up state before tunneling into 
the superconductor. 

We note that the total conductance G — Ga + Gin in 
Eqs. ([l]) - (^) is also an asymmetric function of volt- 
age due to the inelastic processes in the ferromagnet. 
This asymmetry is related to the difference — IIi in 
the maximal cross sections of the majority and minority 
Fermi surfaces and can be used to extract the magnon- 
assisted (inelastic) contribution from transport measure- 
ments. The asymmetric contribution to the conductance 
may be related to the magnon density of states by tak- 
ing the derivative with respect to voltage of the difference 
G{y)-G{-V): 



where F(x) = a;/(l — x^)^. An extreme example occurs in 
tunneling involving a half-metallic ferromagnet which is 
a material where the exchange spin splitting between the 
majority and minority spin bands exceeds the Fermi en- 
ergy, measured from the bottom of the majority band, 
so there, are only majority spin carriers at the Fermi 
energytZI. Since IIj^ = 0, the contributions of conven- 
tional Andreev reflection, Eq. (^, and of states near the 
minority Fermi surface to the inelastic process, second 
term in Eq. (|2|), are absent and the contribution of ma- 
jority states to the inelastic process, first tecai in Eq. (||), 
dominates the subgap /(!/) characteristicsEZl. At the be- 
ginning of Section III we discuss junctions involving a 
half-metallic ferromagnet in more detail and we present 
results (Eq. (41)) for the differential conductance of both 
uniformly transparent and disordered planar interfaces. 

Since the inelastic conductance (||), (^ is written 
in terms of the local density of states ^(uj) of long- 
wavelength magnons at the interface, the same expres- 
sions may be used to describe ferromagnetic systems with 
more complicated magnon spectra. For example, a fer- 
romagnetic superlattice or a trilayer (see, e.g., Ref. ^ 
and Refs therein) can be used as an F element of the 
contact. A suggested geometry, sketched in Fig. 2, con- 
sists of two superconducting reservoirs which sandwich 
a ferromagnetic-normal metal multilayer adjacent to a 
'long' normal layer. The 'long' normal layer produces 
asymmetry in the structure (e.g. by introducing an addi- 
tional resistance) so that the inelastic contribution to the 
conductance would depend on the local density of states 
of magnons at the interface between the multilayer and 
the superconducting reservoir on the left. Note that the 
use of superconducting reservoirs is common when mea- 
suring the current perpendicular to plane (CPP) giant 
magnetoresistance of multilayer structures because they 
reduce the overall resistance of the system. In multilayer 
structures the specteum of magnons may consist of sev- 
eral separate bandsEj so in contacts where the supercon- 
ducting gap energy A is larger than the typical magnon 
bandwidth the magnon band structure may manifest it- 
self in the form of the /(y) characteristics. The precise 
form and the influence of the magnon spectrum in the for- 
mation of the subgap current should depend on whether 
the coupling between magnetic layers is ferromagnetic or 
antiferromagnetic. 

The paper is organized as follows. In Section || we 
introduce the model and technique used for describing a 
tunnel FS contact and give general results for the subgap 
current and the l^V) characteristics valid for an arbitrary 
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interfacial quality. In Section HI we calculate the differ- 



ential conductance for both a half- metallic ferromagnetic 
electrode and a conventional ferromagnet. The results 
for uniform and disordered interfaces are discussed. Ap- 
pendix ^ is devoted to detailed calculations of the subgap 
current based on the Kcldysh formalism. 




F N 




FIG. 2 - Sketch of a suggested geometry in which the 
contribution of magnon-assisted Andreev reflection to the 
differential conductance is related to the magnon density 
of states of a magnetic multilayer system. Two supercon- 
ducting reservoirs (S) sandwich a ferromagnetic-normal 
metal multilayer (F/N) adjacent to a 'long' normal layer 
(N). 



conduction are caused by different groups of electrons, is 
employed to describe the ferromagnetic electrode. The 
magnetism originates from the inner atomic shells (e.g., 
d or f) which have unoccupied electronic orbitals and, 
therefore, possess magnetic moments. On the other 
hand, the conduction is related to electrons with dclocal- 
ized wave functions. These two electron groups are con- 
sidered to be coupled via an intra-atomic exchange intefe 
action. Using the Holstein-Primakoff transformatiorH 
the operators of the localized moments in the interaction 
Hamiltonian can be expressed via magnon creation and 
annihilation operators At low temperatures, where 
the average number of magnons is small < ll^h ><C 25* {S 
is the spin of the localized moments) , the Hamiltonian of 
the ferromagnet Hp can be written as follows 



Hp + Hr, 



Hp 



(7) 



ife = ^ ea(k)4„afe„, e„(k) = e;^(k) - «r7(k)/2, 
a = ± or r,i, (8) 



koL 



(9) 



II. SUBGAP TRANSPORT THROUGH A 
FERROMAGNET-SUPERCONDUCTOR 
INTERFACE 

A. Model 

In what follows we consider a tunnel contact of area A 
between a ferromagnet and a superconductor which have 
dimensions Lp and Ls, respectively, in the direction z 
perpendicular to the contact plane. The contact is posi- 
tioned at z = so that —Ls < z < Lp. The difference 
Es—Ep = eV in the chemical potentials of the electrodes 
is assumed to be smaller than the superconducting gap 
energy A: \eV\ < A and the temperature is finite, but 
also less than the gap energy: T < A. 

To describe electronic tunneling through the barrier at 
the interface {f^ 0) we make use of the tunneling Hamil- 
tonian metjjipcr^ and the nonequilibrium Green functions 
techniqucLjEj. In this approach the total Hamiltonian 
of the system is 



H — Hs + Hp + Hp, 
Ht = Efcfc'a tsF{^', k)4,„afc„ + h.c. 



(5) 
(6) 



where Hp is the tunneling Hamiltonian, Hp{a\^,a}^a) 
is the Hamiltonian of the ferromagnetic electrode, and 
Hs{c^k'a^^k'a) is the Bogolubov de Gennes Hamiltonian 
of a superconductor (hereafter, electron wave vectors in 
the S sidfi-are marked with prime). The so called s-f (s-d) 
modela'E^, which assumes that magnetism and electrical 



Hem = - ,2SY/2 \y<ikpK"^k]<^Pl + 

qkp 

Ugkp = Tu{k,p){ALpy/' J drx;{r)n{r)%{r). (10) 

The first term Hp in (|^) deals with conduction band elec- 
trons which are split into majority e|(k) and minority 
e|(k) subbands due to the s-f (s-d) exchange, ep(k) is 
the bare electron energy, F sets the scale for the splitting 
energy and the dimensionless function 7(k) takes into 
account a possible anisotropy of the band splitting. The 
Hamiltonian iJ.,„ (^ describes free magnons with spec- 
trum LUq which in the general case has a gap ujq=o = too- 
The third term Hpm in (|^) is the electron- magnon cou- 
pling (^^ resulting from the intra-atomic exchange in- 
teraction between the spins of the conduction electrons 
and the localized moments. The matrix element of this 
interaction Uqkp is, generally speaking, of the same or- 
der as the band splitting. In the expression for Uqkp (see 
Eq. (|l^)) <&fe(r) and x<?(r) are the electron and magnon 
wave functions, respectively, which are normalized by the 
volume ALp of the ferromagnet. The dimensionless func- 
tion it(k, p) introduces an extra dependence on the elec- 
tron wave vectors which may come from an anisotropy of 
the electron-magnon coupling near the Fermi surfaces of 
the majority and the minority electrons. 

The tunneling Hamiltonian Hp (|^), which couples the 
two electrodes, is written in terms of the creation and 
annihilation Fermi operators , c and , a of the S and 
F electrodes respectively. The tunneling matrix elements 
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tsF(k',k) describe the transfer of an electron with wave 
vector k from the ferromagnet (F) to the state with k' 
in the superconductor (S). We will consider tsrO^' to 
be a symmetric matrix of the form 



tsF{k',k) = ik'||,k| 



/^2^;|(k; fc;>|,(k||,fc. 



2LsLp 



1/2 



(11) 

(12) 



where 



This expression explicitly takes into account the fact that 
the quantum flows Vg p/ Lg^p perpendicular to the inter- 
face are non-zero on either side of the junction for elec- 
trons participating in the tunneling process; vg^p are the 
perpendicular velocities and es,i?(k) denote the electron 
spectra in the superconductor and the ferromagnet. The 
wavenumbers k are separated into components perpen- 
dicular kz and parallel k|| to the interface, so that the di- 
mensionless factor tk'||,k|| depends only on the momenta 
parallel to the interface. We also suppose that spin is 
conserved upon electronic transfer across the interface 
and pay attention to the spin index a — {t, J.} (or ±) in 
all calculations. 

In what follows it is assumed that the voltage drop 
V occurs across the junction so that both electrodes are 
almost in equilibrium. The total current I through the 
interface is proportional to the average value of the rate 
of change of the number of particles in, for example, the 
ferromagnet: I = e < Np >. This rate can be found 
from the commutator of Np = X^fca '^la'^fea with the to- 
tal Hamiltonian H = Hg + Hp + Ht- Bearing in mind 
that only the tunneling Hamiltonian Ht fails to commute 
with Np, we arrive at the following formula for / 



kk' a 

't*gp{k',k)< 



(13) 



The averaged quantities in Eq. (y_3|) can be expressed 
in terms of noncquilibrium. (Keldysh) Green functions 
g+-{t,t) and 0-+(t,i)El"El and it is convenient for the 
case of a superconducting electrode to use the 2x2 
Nambu representation as shown explicitly in Eq. (|Al|) of 



Appendix The anti-hermitian relationships Eq. (A2) 
for and Q ^ enable us to represent the current / 
( p^ in terms of two off-diagonal (with respect to the 
electrode subscripts {F, S}) Keldysh functions: 



2eReY, [tsF(k', k)g+-^^(kt, k't)- 

kk' 

+t*gp{k','k)gp+^^^{kt,'k't) , 



(14) 



where the subscripts {1,2} refer to spin indices and the 
symbol Re means the real part of the sum. 



B. Andreev and magnon-assisted currents. 
Derivation of the subgap I(V) characteristics 



The Green functions Qpg^-^^ and Gpg^22 ( P^ can be 
determined with the help of perturbation theory. For 
this purpose, we represent the total Hamiltonian (^) in 
the form H = Hs + He+ Hm + {Hp + Hem) where the sum 
of the tunneling Hamiltonian and the Hamiltonian of the 
electron-magnon interaction is treated as a perturbation. 
In this approach the usual Andreev corrections to Green 
functions Gpg^n and Gpg^22 appcar ui the third order of 
the perturbation with respect to the tunneling (and ze- 
roth order in the electron-magnon interaction) and are 
proportional to 



/ dTldT2dT3HT{Tl)HT{T2)HT{T3), 
J K 



(15) 



where the integration is performed along the Keldysh 
contour (see Eq. (A3) in Appendix ^). They give rise 
to a tunneling current / ( |T^ which is of the fourth order 
in the hopping elements tk' ^k as it should be because two 
particles are involved in the process of Andreev reflec- 
tion. The lowest order contribution to magnon-assisted 
Andreev reflection also arises in the third order of the 
tunneling matrix tk'^ki but it is bilinear in the electron- 
magnon interaction. 



/jj, dTldT2dT^dTidT^ X 
X Ht{ti)Ht{t2)Ht (ts ) 7?em {TijHem (rg ) . 



(16) 



Since the interaction Hamiltonian H^m (^0|) is of the first 
order in the magnon operators 6^^ bq, the diagrams cor- 
responding to this term will contain only one magnon 
line. The diagrammatic calculation of the tunneling cur- 
rent ( p^ ) performed in Appendix ^ show that / can be 
represented as the sum of three different contributions: 



'TT 



(17) 



4e 



'Ti 



fel,fe2 



^ tgp{]^M)tsF{~\^. k2)gss|i2(k'; 



5{~e + eV- £.kiA)5{2eV - ^ij - ^2,1) x (18) 
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q,kl,k2 



J2 isF(k', ki)isj.(-k', k2)5ss|i2(k'; 



E 



Uqk2p 



(19) 



X(5(tJ - LUq)d{-e + eV- ^ki,])Si2eV ~ ikl^ ~ Cfc2,T - ^) {"(-Cfci,T)"("62.T)[l + N{uj)] - n{^kiA)n{^k2,]')N{uj)} , 



J2 tsFik', ki)tsFi-k', k2)5ss|i2(k') 



E 



U, 



qk2p 



(20) 



x,^(a; - Ljq)5{e + eV~ Ckui)Si2eV - ^ki,i - ^k2,i + ^) {ni-^ki,i)n{-ik2A)N{ij) - n{^ki,iM^k2,i)[l + N{ij)]} . 



The first term in Eq. ( |l7|) is the usual Andreev current 
III which, as can be seen from Eq. (|l8|), results in ei- 
ther the gain or loss in the ferromagnet of two carriers 
belonging to different spin bands. Their energies S.ki.i 
and ^fe2,i measured with respect to the chemical poten- 
tial on the F side Ep obey a conservation law of the 
form ^ki,i + ^fc2,i = 2eV where 2eV is the energy of 
a Cooper pair with respect to Ep. The current In is 
non-zero owing to the imbalance between the gain of the 
two electrons, n{-^ki,f )n{-^k2d) = [1 - "(Cfci,T)][l 



"•(Cfc2,x)], and their loss, n{£^ki,i)n'{^k2,i)- In equilibrium 
{V = 0) energy conservation ensures vanishing I^. The 
factor I J2k' *SF(k',ki)ts_F(-k',k2)5ss|i2(k')p describes 
the two-particle tunneling process resulting in the cre- 
ation or annihilation of a spinless Cooper pair with an 
anomalous average gss|i2(k') on the S side. Using the se- 
lection rules expressed by the delta- functions in Eq. (|T^), 
one can find that at zero temperature, T — 0, the An- 
dreev current (|l^) assumes the following form 



J 



Sire 



eV 



de 



eV 



tsFik', ki)tsFi-k', k2).9s5|i2(k', e) 



(5(-e + eV- ^kiA)S{2eV - Cfci,T " 62a) 



(21) 



The magnon-assisted contributions to the subgap cur- 
rent ( |l7|) are given by Eqs. ( [l^ ) and (^. Let us an- 
alyze first the contribution /|| from the majority carri- 
ers. In Eq. (19) the product of the occupation numbers 
"-(~Cfci,T)'^(~^fc2,T)][l + represents the gain of two 

spin up electrons and a magnon in the F side. This pro- 
cess can be viewed as the injection of two electrons with 
energies ^ki,i and ^p^i into the ferromagnet where the 
spin down electron has a finite lifetime because of the 
intra-atomic exchange interaction with the localized mo- 
ments. Then the spin down electron emits a magnon of 
energy uj and occupies an empty state ^fc2,t in the spin 
up conduction band. Note that the energy in the de- 
nominator of Eq. (p^), ^p_| — £^k2,i — is related to the 
inverse lifetime of the spin down electron in the virtual 
state. The reverse kinetic process consists of the loss of 
two spin up electrons and a magnon from the F region: 
'>T-{^ki,'\)n-{^k2,'\)N{uj). It is possible only at finite temper- 
atures because one of the spin up carriers has to absorb 
a magnon in order to form the necessary intermediate 
spin down state before tunneling into the superconduc- 
tor. Both processes obey an energy conservation law of 
the form 2eV = Cfei.T +?fe2,T Thus, if the energy 2eV 
of a Cooper pair is equal to zero (at V = 0), the direct 
and reverse processes compensate each other. 

At zero temperature, where the majority current In 
( |l9| ) is determined by the gain of electrons in the F side, 
the integration over the electron energy e is reduced to 
the interval uj — eV < e < eV. As a consequence, only 



magnons with energies less than the energy of a Cooper 
pair, i.e. < a; < 2eV, contribute to the current. Fur- 
thermore, the last inequality shows that the majority 
current is non-zero only at positive voltages. The min- 
imal value of V is related to the gap energy ujq in the 
magnon spectrum: LUQ/2e < V. On the other hand, the 
inequalities uj < 2eV < 2A imply that the magnons of 
interest are mainly long-wavelength ones because the su- 
perconducting gap energy A is usually much less than the 
magnon Debye energy. In other words, in Eq. ( |l0|) for the 
matrix element Uqk2p of the electron-magnon interaction 
the magnon wave function varies much more slowly than 
the electron ones. For this reason, Uqk2p falls rapidly as 
the difference |p — k2| in the electron momenta increases, 
which means that both the virtual (spin down) and the 
final (spin up) electron states are close to the majority 
Fermi surface. It allows us to neglect the change in the 
electron momentum in the denominator of Eq. ( p^ and 
calculate the sum over the momentum of the virtual state 
p as follows 



E 



Uqk2p 



1 



- Cfc2,T - ^ r7(k2) 



(22) 



-7(k2)^ 

In order to obtain the last equation we have used for- 
mula ( pO| ) for Uqk2p with wave functions for electrons in 
a large volume with a hard wall at z = 0, i.e. ^^(r) = 
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2^^'^{ALp) exp(zk||r||) sin(fczz), and magnon wave 
functions of the form Xq(r) = A~^^'^ exp{iq\\r\\)xq^{z), 
where Xq^ (^) corresponds to motion in the direction per- 
pendicular to the contact plane. It can be shown that the 
integration over p in Eq. gives rise to a <5-shaped fac- 
tor which in the limit (kz/qz) oo transforms into S{z). 



That is why the effectiveness of the electron-magnon in- 
teraction depends on the magnitude of the magnon wave 
function near the interface (Eq. (p2|)). According to (|2^), 
the contribution of Eq. (|l9| ) can be represented for T ^ 
as 



J 



Jo Juj-eV 1^^1.2 



^ tSF(k', ki)isF(-k', k2)555|12(k') 



(23) 



where Vl{u}) denotes the magnon local density of states 
at the interface 



(24) 



Since the magnon wave function Xq^ {^) normalized 
on the length Lp of the ferromagnet, the quantity 
Lp\xq^ (O)l^ in Eq. ( |2^ ) is the dimensionless magnon am- 
plitude at the interface which does not depend on Lp. 

The other magnon-assisted contribution ( pO| ) is 
formed by minority electron states. The kinetic processes 
giving rise to consist of the gain or loss in the F re- 
gion of two spin down electrons. In the direct process (the 
gain) Andreev reflection results in the injection of both a 
spin up and a spin down electron with energies and 
^fei,i respectively. Then the former absorbs a magnon 



of energy uj and occupies a state £,k2,i in the minority 
conduction band. This takes place only at finite temper- 
atures. In the reverse process (the loss) one of the minor- 
ity electrons emits a magnon and forms an intermediate 
spin up state with energy This allows it to tunnel 
into the superconductor along with another spin down 
carrier having energy ^ki.i- Both processes obey the en- 
ergy conservation law of the form 2eV = S,ki,\, +^fe2 j, — 
It can be shown in an analogous way as for J|| ( psj ) that 
only interactions with long-wavelength magnons are of 
importance, which take place in the vicinity of the mi- 
nority Fermi surface. We note that, in contrast to the 
direct process, the reverse one is possible at zero temper- 
ature. That is why the T — > asymptotic form for the 
current ( po[ ) is determined by the 'loss' term in Eq. ( |20| ) 
and, therefore, it is non-zero only at negative voltages 
V < -uJo/2e: 



J 



hi = -T^0(-e^ - ^o) / dujniuj) / de V 

Jo Ju^+eV ;,,fe2 



^(k2)P 
7^(k2) 



^tsF(k',ki)tsF(-k',k2) .9SS 1 12 ) 



(25) 



In Eqs. (^), ( p3| ) and ( psf ) the anomalous Green func- 
tion 5(ss|]^k', e) oLthe superconductor is given by the 
well-known formulaeJ (see also Eq. (A19)) which can be 
rewritten as follows 



SSS|i2(k',e) 



A 



(A' 



2 _ ,2Nl/2 



(A2 _ £2)1/2 _ Es)^ + A2 - £2 



(26) 



Here the second multiplier represents the Lorentz curve 
as a function of e5(k') — Es for subgap excitations with 
e < A. In the Andreev approximation, where A ^ Eg, 
it takes the delta-shaped form 



5ss|i2(k',e) 



ttA 



(A2 



— S{Es-es{k')), 



(27) 



that allows us to decouple the integrations over the en- 
ergy e and the electron wave vector k' in Eqs. (p|), (12 



and (|25|). This means that the I{V) characteristics of 
the contact turn out to be independent of the detailed 
structure of the interface. Also, taking into account that 
e| < \eV\ -C r^Es^Ep, we have for the usual Andreev 
current 



4e ^ , eF 

7ti = — Aarctanh— — x 
h A 



(28) 



J rfP||l y C?P||2-Ra(P||1,P||2), 



(29) 



Here H| and H| are the areas of the maximal cross sec- 
tions of the majority and minority Fermi surfaces in the 
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plane parallel to the interface. The values p-\^max and 
Pl,max of the perpendicular momentum correspond to the 
positions of the maximal cross sections in the p-space. 
The matrix i?A(P||i, P||2) in ( ^8|) describes Andreev re- 
flection of a spin up channel specified by the value of 
parallel momentum p||i into a spin down channel with 
parallel momentum p||2 or vice versa. The integrals in 
Eq. ( p8| ) take into account all possible transitions between 
the spin up and spin down modes in the contact. The 
form of the reflection matrix depends on the quality of 
the interface and is given by 

ra{p\\i,p\\2) = (^-^^ I J dpii t(p^|,P||l)^(-p^|,P||2)l^ 



(30) 



)), (31) 



where Us is the area of the maximal cross section of the 
Fermi surface in the superconductor. The expression for 
the reflection matrix ( ^0| ) can be derived by integrating 
over the perpendicular momenta in Eq. (21) with the use 
of the delta- functions and the formula (11) for the tun- 
neling matrix elements. 

The inelastic contributions to the subgap I{V) char- 
acteristics are given within the Andreev approximation 
by 



'TT 



A 

hS Kh^ 



J dp\\l y dP||2^m(P||l,P||2) X 



(32) 



X Ae{eV - Wo) / n{u;)duj 
Jo 



2eV 



eV eV -u 

arctanh— — h arctanh : 

A A 



J =_(A 
hS \ 



J dP\\l J dp\\2Rm{P\\l,P\\2) ^ 



(33) 



X Ae{-eV - ujq) / n{uj)duj 



-2eV 



eV eV + uj' 

arctanh— — h arctanh : 

A A 



-Rm(P|ll,P||2) 



l^(P||2)P 
7^(P||2) 



j fiPii *(p'|,Plli)i(-p'|,P|| 



ns 



(34) 



Since a spin-flip process associated with magnon emis- 
sion compensates the change in spin resulting from An- 
dreev reflection, the matrix i?in(P||i, P||2) in Eqs. (^) 
and (^3|) for the magnon-assisted currents describes 
transitions between electron channels of the same spin. 
Note that i?m(P||i, P||2) is proportional to the factor 
■"(P||2)P/7^(P||2) that takes into account an anisotropy 
of the intra-atomic exchange interaction near the Fermi 
surface of the majority and minority electrons. 



C. Uniformly transparent interface 



Let us analyze now the properties of the matrices (^) 
and (^) for a plane, wide contact with a clean barrier 
where electron momenta parallel to the interface are con- 
served upon tunnehng. In this case the factor i(p'|,p||) 
in Eqs. ( pO| ) and (jsj) is proportional to the matrix ele- 
ment of a transition between two plane waves normalized 
by the contact area A: 



i(p'|,P||) = t 



/iW-pil) 



A 



(35) 



where \t\^ represents the tunneling probability per chan- 
nel. We assume that |ip is much less than unity and is 
related to the barrier width which is, in turn, considered 
to be the same for all points along the contact. For a con- 
tact with large, but finite area the peak of the function 
'^(p'l ~ P||) at p'l = p|| is as high as the ratio A/h? . For 
this reason the factor /i^(5(p|| —p\\)/A in Eq. ( |35| ) should 
be treated as a Kronecker delta. Assuming for simplicity 
that the maximal cross section Us of the Fermi surface in 
the superconductor to be the largest one: 115 > > Ifj, 
we obtain matrices (30) and (p3) in the form 



^/i2j(p||i +P112) 



^a(P||1,P||2) = 1^1' ^ 
^Jn(P||l,P||2 = Ra{P\\1,P\\2)^-, ^• 

7 (P1I2) 



(36) 



As it should be for a homogeneous interface, the rela- 
tionship between the parallel momenta of the incident 
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and Andreev reflected particles on the F side is the same 
as for electrons forming a Cooper pair on the S side: 

P||l +P|12 = 0. 

D. Randomly transparent interface 

Now we introduce a model for describing a strongly 
nonuniform interface which is transparent only in a fi- 
nite number of points randomly distributed over an area 

— 1/2 

A with density n^. The typical distance between 
the transparent points is assumed to be much larger than 
the Fermi wavelengths in the superconductor A5 and the 
ferromagnet A^: 



As,F"l^' « 1. 



(37) 



Each transparent point contact can be treated as a de- 
fect which causes electron scattering in the plane parallel 
to the interface. Under condition (|37| ) we can describe 
the system of such defects by a sum of short-range po- 
tentials a,J2j ~ where Vj is the position vector 
of the jth contact in the plane, a A5 ^ is its area, and 
tj introduces the magnitude of the potential at r = Vj. 
The transfer matrix element i(pj| , p||) (see Eqs. ( |30| ) and 
(p^)), which accounts for the change in the parallel mo- 
mentum of an electron upon tunneling, should be identi- 
fied as the matrix element of the total scattering potential 
calculated with the use of plane waves normalized on A: 

^(p'i,Pii) = ^^Xl*j'^('"~''j)'^^p[*^"^(P|| ^pIi)*"] 



shift is equal to zero. As a result, the contribution from 
such terms is (a/A)'^^^ l^jl^- Next, one has to average 
the product of the phase factors in (^) divided into all 
possible pairs rji — rj2,rj3 = rj4, ... It is easy to show 
that the contribution from these terms is small compared 
to the contribution mentioned above by the parameter 
{Xs-F^nh -C 1. Averaging of the four independent phase 
factors in Eq. ( p9| ) gives an even smaller result. Thus, the 
averaged reflection matrices (|30|) and (^ih are 



(40) 



i?4P„.P„2) = (^)Y^)'El^.l 



-Rm(P||i,P||2) = Ra 



\ J 

^Kpi|2)P ^ 



7^(P|12) 



, i?A < 1. 



(38) 



The matrix Ra for the usual Andreev process is inde- 
pendent of the momenta and proportional to the squared 
number of transverse electron modes alis/h^ '--^ 1 in one 
point contact on the S side. The factor |tj|^ gives 
the total transparency of the interface summed over all 
the contacts. Since we have treated tunneling as a per- 
turbation, each tunneling probability \tj\'^ must be much 
less than unity. Note that the momentum dependence of 
the inelastic scattering matrix i?i„ is related only to the 
anisotropy of the electron-magnon interaction. 



III. DISCUSSION OF THE MAGNON-ASSISTED 
CONDUCTANCE OF A 
FERROMAGNET-SUPERCONDUCTOR 
JUNCTION 



Here the constants tj are chosen to be dimensionless. 



Since the reflection matrices (|30| ) and (|3^) now become 
random functions of Vj , we have to average them over the 
position of each defect in the same spirit as in the prob- 
lem of electron-impurity scattering in bulk metals (see, 
e.g., Ref. ^ and Refs therein). The product of matrix 
elements to be averaged (see Eqs. (^ and (p^)) is 



i(p'|l,P||l)^(-p'|l,P||2)i*(P||2iP||l)^*(-P||2'P||2) 
I* 
A^ 



ilj2j3j4 



■J4 



X 



xexpi?i ^ (p||i - p||j)rji + (p||2 + p||Jrj2- 



-(P||l -P||2)rj3 - (P||2 +P'|2)rj4 



(39) 



where the bar represents averaging. The sum in the argu- 
ment of the exponential function is the total phase shift 
caused by scattering at the contacts located at points 
rji,rj2,rj3 and rj4. If these points correspond to the 
same defect, i.e. Vji ~ Yj2 = r^a = rj4, the total phase 



The band structure of a ferromagnetic electrode near 
the Fermi level may have a dramatic effect on the sub- 
gap conductance of an FS contact. Let us demonstrate 
this by considering the extreme example of a junction 
involving a half-metallic ferromagnet. 

In half-metallic ferromagnets the splitting F between 
the majority and minority subbands is so strong that 
the bottom energy r7(0)/2 of the minority subband (see 
Eq. (||)) is higher than the Fermi level, i.e. r7(0) > 2Ef. 
For this reason only majority carriers are present at the 
Fermi level, and, therefore, usual Andreev reflection is 
forbidden by the s pin current conservation law. Indeed, 
according to Eq. (p8[), the usual Andreev current van- 
ishes as the number of the minority spin channels Alii/h^ 
goes to zero. The same is true for ( |3^ ) because this 
magnon-assisted contribution is determined by the mi- 
nority states near the Fermi energy. Thus, the current 
in the half-metallic electrode is completely spin polar- 
ized (mediated by the spin up carriers) and its value / is 
equal to the inelastic magnon-assisted contribution (IsT 
/ = Iw . 
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FIG. 3 - Tunneling process between a superconduct- 
ing electrode on the left hand side and a half- metallic 
ferromagnet on the right for {a) V > Q and (b) < 0. 
For (a) F > a down spin electron tunnelling into the 
fcrromagnet may emit a magnon and incorporate itself 
into the majority (up) conduction band. For ih) V < Q 
no spin flip process is possible at T = because in the 
initial state (right) there are no thermal magnons for an 
up spin electron to absorb. 



dpil |u(p 



7^(P|l) 



According to Eq. (|32|), the liV) characteristics of the 
contact are asymmetric with respect to the bias voltage 
V . The current is zero for negative bias and finite for 
V > wo/2e. This asymmetry can be explained using 
the sketch in Fig. 3 which shows the tunneling process 
between a superconducting electrode and a half-metallic 
ferromagnet for (a.) V > and (b) V < 0. As discussed in 
the Introduction, for V > 0, Fig. 3(a), Andreev reflection 
manifests itself as the injection of both a majority and a 
minority electron into the ferromagnet where they have 
energies eV±e with respect to the Fermi level Ep- Since 
there are no spin down states near the Fermi level in a 
half-metallic ferromagnet, the state of the minority elec- 
tron turns out to be virtual and the process of magnon 
emission allows this electron to occupy an empty state in 
the majority conduction band above Ep. 

For V < 0, Fig. 3 (b), Andreev reflection would re- 
sult in the injection of a singlet pair of electrons from 
the feromagnet into the superconductor. With respect 
to the Fermi level in the ferromagnet Ep these energies 
are — |el^l -|- e and — leVj — e, respectively. A spin up 
electron cannot emit a magnon due to conservation of 
total spin in the exchange interaction so that the only 
possibility could be that a spin up electron in the ferro- 
magnet would absorb a magnon to form an intermediate 
spin down state before tunneling into the superconduc- 
tor. However there are no thermally excited magnons 
at T = in the initial state of the ferromagnet, so it 
is impossible for magnon-assisted Andreev reflection to 
contribute to current formation in negatively biased con- 
tacts. 

From Eqs. (||), (|^) and @) we flnd that the 
inelastic differential conductance Gm of a half-metal- 
superconductor contact is 



E, It. 



\t\^iAU,/h^) 
^{aUs/hy{aU^/hY 



ideal interface 
disordered interface 



(41) 



where the voltage-dependent factor M{x) caused by 
the inelastic processes in the ferromagnet is defined in 
Eq. (H). Note that in the case of a disordered interface, 
the conductance ( ^ ) is proportional to the product of 
four electron channels (two on either side of the con- 
tact) participating in transport through one point con- 
tact. This means that, in the absence of the parallel 
momenta conservation law, two electron channels from 
one side of the junction are independently transformed 
into two channels on the other side. 

In the case of conventional ferromagnets the subgap 
current / contains all three contributions (21), (E3h and 



3 



4 (aTls\ ( aHt 



M{eV) 



(p5|), and, therefore, the differential conductance has 
both the elastic (Andreev) and inelastic parts: G — 
Ga + Gin- For tlie case of a uniform interface, Ga and 
Gin can be obtained using Eqs. (^), (p3|), ( |25| ) and Eqs. 
( p6| ) for the reflection matrices Ra and Rin- The results 
are given by Eqs. (|l|) ~ (^. Calculating Andreev and 
magnon-assisted conductances for a disordered interface, 
we find that 



Ga{V) 



4e^ 



E 



dV\\ l"(P||)l' 
Ht 7^(P||) 



M{-eV)^ 



M\ l"(pii)l' 



HtJ Ht 7'(P||) 



(42) 



(43) 
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In contrast to the case of a half-metalhc ferromagnet (see 
Eq. in the general situation the inelastic conduc- 

tance^^) is non-zero at negative bias as well. As dis- 
cussed in the Introduction and Section 0, this feature is 
explained by the presence of magnon-assisted transport 
of the minority carriers. 

We note that the electron and magnon 'contributions' 
to the magnon-assisted conductance Gin fall into sepa- 
rate multiphers (see, e.g., Eqs. (|4l| ) and (|4^)). On the 
one hand, this fact is due to the small momentum trans- 
fer of the electron-magnon interaction (see Eq. (p^)). On 
the other hand, it is necessary that the superconducting 
gap energy is much less than Fermi energy, A <^ Eg , (see 
Eq. ( p7[ ) for the anomalous average under the Andreev 
approximation). Within these reasonable restrictions our 
approach allows us to calculate the tunnel conductance 
for arbitrary electron and magnon dispersion spectra and 
an arbitrary shaped potential barrier at the interface. 
Although the interfacial disorder does not influence the 
form of the I{V) characteristics (see Eqs. (^T|), ( ^ and 
(p5|)), it determines the dependence of the conductances 
on the parameters of the Fermi surfaces in the super- 
conductor and the ferromagnet (cf. Eqs. (^, (|^) and 
Eqs. (H), (H)). 

According to our analysis, the inelastic processes of 
magnon emission should manifest themselves in the 
asymmetric behaviour of the differential conductance 
G{V), which is a direct consequence of the spin polariza- 
tion of carriers in the ferromagnetic electrode. In order 
to observe the inelastic contribution to G{V) it is neces- 
sary to eliminate the conventional Andreev (symmetric) 
part of the conductance (|l|), (^2|). Besides, the proposed 
theory demands that the gap in the magnon spectrum 
must be smaller than the superconducting gap energy, 
< A. The latter requirement might not be satisfied in 
the experiments of Refs. [l0|-p^. That is why the results of 
these works can be treated neither in favour nor against 
the existence of the magnon-assisted Andreev processes. 
On the other hand, the requirement < A can be met in 
soft half-metallic materials, e.g., manganiteal^. In Refs. 
p|,p^ these compounds were used in tunnel contacts with 
high-Tc superconductors (materials with large A). The 
(dl/dV) — V curves measured in Ref. ^ clearly demon- 
strated asymmetry for |V^L < A/e which is consistent 
with that shown by Eqs. (g) and (^). As pointed out 
in Ref. the asymmetry of the (dl/dV) — V curves 
might be attributed to the d-wave type pairing in the 
high-Tc superconductor. However d-wave simulations 
performed in Ref. |l^ gave good results for quite large 
voltages \V\ « A/e, but did not explain the asymme- 
try in the range \V\ < A/e where the magnon-assisted 
Andreev reflection may occur. 
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N. R. Cooper, D. M. Edwards, C. J. Lambert and 
Yu. V. Nazarov for discussions. This work was supported 
by EPSRC, NATO and the Royal Society. 



APPENDIX A: KELDYSH FUNCTIONS 
FORMALISM FOR CALCULATING SUBGAP 
TUNNELING CURRENT 

Following Ref. we introduce the 2x2 Nambu repre- 
sentation for the nonequilibrium (Keldysh) Green func- 
tion in the system of units where h — 1: 

^'A'M^ (Al) 

V < TKC^i(^j)4T(^J') > < TKcl^{t,)cE'i{t'^,) > ) ■ 

Here Tk is an operator ordering the times along the 
Keldysh contour which consists of a "positive" branch 
going from — oo to -|-oo and a "negative" one going from 
-f oo to —oo. Any time of the negative branch is consid- 
ered as posterior to any time of the positive branch. The 
index J (j') means "-I-" or "— " depending on which branch 
of the Keldysh contour the time argument tj {t'j, ) belongs 
to. The subscripts E and E' label the electrodes (S or 
F) in which a particle is created or annihiliated. Using 
these definitions, it is easy to check that, for instance, the 
Keldysh function Q^^^.^^ = — i < Txat(t)a|(t') > con- 
tains chronological ordering of the operators and, thus, 
it is the usual causal Green function of the spin up elec- 
trons in the ferromagnet. We also note that according to 
Eq. (^ the functions e^F>a'(*'*') and Qst'laa'^^^^') 
are anti-hermitian, i.e.. 



tei«.'(*.^')* = -e+rF|„,„(t',0, 



(A2) 



Integration over the Keldysh contour is reduced to the 
interval —oo to -\-oo as follows 



/ dt...= 

Jk J-oo 



dt J2 ^■ 

j= -1-1,-1 



(A3) 



Since we deal with a static problem, it is convenient 
to go over to the energetic representation for Green func- 
tions Q^g]^ii and ^^g|22 in (0)- Then they can be calcu- 
lated using conventional perturbative diagram technique 
(see, e.g., Ref. ^) where the perturbation Hamiltonians 
are given by Eqs. (^) and ([iq). If we multiply diagrams 
for Q'^'^^ii and Qp^\22 ^spi}^ ^ k) and tj^(k', k) respec- 
tively and integrate over all momenta and energies, we 
will obtain diagrams in the form of closed loops which 
represent the perturbation series for the tunneling cur- 
rent / (|lj). These diagrams are shown in Fig. 4. 

Solid lines in Fig. 4 represent unperturbed electron 
Green functions of the ferromagnet 5^^111' 5ff|22 and 
the superconductor 5'55|]^2' 555121 ■ flashed lines are the 
unperturbed magnon functions d^^ . The arrows indicate 
the order of the arguments (momenta and Keldysh in- 
dices = +,— ). Each tunneling matrix element tsF 
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couples Green functions of the different electrodes. In 
other words, two solid lines connected by tsF represent 
an elementary hybrid block which gives one Nambu ma- 
trix element of FS or SF Green function to the first or- 
der in the tunneling constant. Since both the FS and SF 
functions arc 2x2 matrices, there are 8 different hybrid 
blocks. Note that the diagrams shown in Fig. 4 contain 
only the off-diagonal in Nambu indices blocks which are 
relevant to the process of Andreev reflection. The elec- 
tron energy e is assumed to be measured from the Fermi 
level Es in the superconductor. Thus electron Green 
functions of the fcrromagnet have e + aeV as the energy 
argument. 



a) 




b) t.,.- 



FFIll . - FI32: 



^ FF|22 ^ FF|I' ^ ^ FFP^i .. 



tsF 




FIG. 4. - Diagrammatic representation of (a) Andreev 
current I^i, (b) magnon-assisted contribution from the 
majority carriers /|| and (c) magnon-assisted contribu- 
tion from the minority electrons In . 

The diagrams in Fig. 4(a) correspond to the usual An- 
dreev process where the incoming electron from the fcr- 
romagnet and the Andreev reflected one occupy the op- 
posite spin bands. As we will see later, the sign "— " in 
Fig. 4 ensures that only causal Green functions enter the 
expression for the current. The diagrams in Figs. 4(b) 
and 4(c) describe Andreev reflection which is preceded 
or followed by spin-flip processes caused by the emis- 
sion or absorption of magnons. These processes compen- 
sate the change in spin resulting from Andreev reflection. 
That is why both the incoming and the outgoing elec- 
trons belong to the same spin band. In the diagrams in 
Fig. 4(b) the initial and final electron states are both spin 



up ones. They are described by Green functions with the 
subscripts FF\11. Green functions _Fi^|22 correspond to 
intermediate (virtual) states of a spin down electron. As 
pointed out in the previous sections, such a state can be 
formed as a result of Andreev reflection, then the spin 
down electron emits a magnon and incorporates itself 
into an empty state in the spin up conduction band. Al- 
ternatively, a virtual spin down state can be formed as a 
result of the absorption of a magnon by a spin up elec- 
tron. Then this spin down state is Andreev reflected 
as a majority carrier. Since the formalism developed in 
this section does not require temperature to be zero, the 
diagrams in Fig. 4(b) take into account both of these 
processes which lead to the existence of 'gain' and 'loss' 
terms in the equation for the majority magnon-assisted 
current (see Eq. ( [T9| ) in Section |l|) . In the diagrams in 
Fig. 4(c) both the incoming and outgoing electrons be- 
long to the minority spin band. They are represented by 
Green functions with subscripts Fi^|22. Such an electron 
can either emit a magnon and form an intermediate spin 
up state, which then experiences Andreev reflection, or 
first be Andreev reflected as a majority carrier which then 
absorbs a magnon. These competing processes result in 
the 'gain'-'loss' structure of the minority magnon-assisted 
current (see Eq. (^0|) in Section |l|). 

Each Green function in Fig. 4 is expressed in terms of 
the retarded g^, advanced and Keldysh functions 
(rf^, d"^ and for magnons) by means of the relation- 
ship 



9' 



(A4) 



which together with ( |A3| ) allows us to formalize summa- 
tions over the Keldysh indices. Indeed, at a given vertex 
labeled by the index j the Keldysh sum has the structure 
, where the factor j™ (m is integer) comes 



+ 1: 



1+m 



from the product of Green functions (A4) connected by 



this vertex. It is clear that only those products of Green 
functions which give odd powers of j 'survive' after the 
summation: 



j=+i -1 



2 for odd m 
for even m 



The Keldysh functions g^ and d^ are related to retarded 
and advanced functions as follows: 

5FF|aa = [9FF\c.ai^ + "eF) " gjpi^Je + acV)] X (A5) 
x[l-2n{e + aeV)], 

d^{u) = [d^iu) - d^iu)][l + 2N{lo)] 

with n and N denoting the electron and magnon occupa- 
tion numbers respectively and uj is the magnon energy. 

After some algebra the tunneling current / (14) can be 
represented as the sum of three different contributions 
(see Eq. (|l7|)) given by 
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[il..i4] 



(A6) X (c'^(q2,qi,w).g^^|22(k2,ki,-e-w-ey)- 
+rf'^(q2, qi, t^)5FF|22(k2, ki, -e - w - ey) + 



9ff\ii^sf9ss\i2^*sf9ff\22^sf9ss\2i^sf 

~9^FF\22^SFg^SS\21^SFg^pp\llt*sFg^SS\12^*SF'^ 



(A7) 



[il..i6] 

is 



21 (* 22 4-* 23 ^^24 25 I 26 J- 

y_FF|ll''SFy55|i2''SFy_FF|22'^22yFF|22''5'-PySS|2l''S'f' ~ 
~ffFF|22^22.9FF|22*'SF355|2i*SF5]^^|ll^SF5ss|i2^SF} 



7£ i: 



(A8) 



[2l..t6] 

J.* U 



5FF|ll^ll3FF|ll'^SF5sS|12%F5FF|22'^5'-P5'sS|21^5'f' ~ 
~9FF\22^SFgss\2V'SFgpp\il^ligpp\lli-SF9ss\12^SF j ; 



Expression (A6) involves the sum over the set of indices 
[h.-iA] = KAAA, RKAA, RRKA, RRRK, i.e. each 
term in this sum contains one Keldysh function muhi- 
plyed by the advanced functions on the right and/or by 
the retarded on es o n the left. The analogous sums in 
Eqs. (|A7| ) and (A8) have six terms of the structure of 
KA..A, RKA..A,.. and R..RK . Such terms have poles in 
both half-planes of the complex variable e and, therefore, 
give nonvanishing contributions to the current. Each 
Green function and the hopping element in the curly 
brackets in Eqs. (A6)-(A8) depends on two momentum 
variables which are considered to be matrix indices (dis- 
crete, in the general case) so that the products in the 
curly brackets are treated as matrix ones. The symbol 
Tr means the sum of all the diagonal elements of these 
products. By T,^^, and we denote the self- 

energies for the electron-magnon interaction: 

= E C^.i,"^Pif^,W/^x (A9) 

9l,92,fcl.fc2 

X (d^(q2,qi,t^)5FF^22(l<^2,ki, -e-uj- eV) + 
+d^(^)(q2, qi, c^)gf;^|22(k2, ki, -e - w - eV)) , 



q2,k2.p2 I ^ X (AlO) 



qi,q2,ki,k2 

X (d^'^(q2,qi,t^)5FF^n(k2,ki, - u + eV)+ 



+d^(^)(q2, qi,^)gfj.|n(k2, ki, -e - cu + eV)) , 



= E U,,,ki..piU;,^k2.^, J ^ X (All) 

9l,l}2,fcl,fc2 



+d (q2, qi, cj)gj,^|22(k2, ki, -e - w - eV) 



qi,q2,ki,k2 

X (^d^{c[2,q.i,uj)gpp^ii{k2,ki, -e - uj + eV)+ 
-Hd^(q2, qi, w).g^p|ii(k2, ki, -e - cj + eV) + 
+rf^(q2, qi, w)5fFlii(k2, ki, -e-uj + eV) 



Sl(e + aeV) - [El(e + aeV) - El(e + aeV)] 
x[l - 2n{e + aeV)]. 



(A13) 



In the subgap regime \eV\ < A Eqs. (|A6D-(|A8|) can 
be simplified. Since there no quasiparticle states in a 
superconductor in the energy range of interest (|e| < 
\eV\ < A), we can set 5fs|i2 = ffssii2 = 3ss|i2 and 
9ss\2i = 9ss\2i = 9ss\2i in Eqs. (0)~(0). In line with 



Eq. ( [A5| ) all terms containing the Keldysh functions of 
the superconductor, ggg^i2 and 3f^|2ij vanish: 



de. 



eRe J —[n{e + eV) - n{e - eV)] X 
xTr|(.g^^|ii ~ gFF\ii)i*SF9ss\ 

x(.9fF|22 ~ 5ff|22)*S-F5SS|21*5f| 



(AM) 



12^SF X 



'TT 



eRe 



de 
2^ 



n(e -I- eV) - n{e - eV)] x 



(A15) 



Tr |(g:J^F|ii ~ gFF\ll)i*SF9SS\12t*SP^ 
X {g^p\ 22^225fF|22 ~ 5ff| 22^22.9ff| 22)^'SF5sS|2lisF | , 

^[n{e + eV)~n{e-eV)]x (AI6) 

|(ffFF|ll^ll3FF|ll ~9FF\ll'^ll9FF\ll)^*SF9SS\12t*SF^ 
X (.9ff 22 gFF\22)^SFgsS\2ltSF^ ■ 

Further simplifications can be made in the case of the 
homogeneous electrodes which are described by Green 
functions of the form (see, e.g., Ref. Eoh 



5FF|oJki,k2,e) - [, + «(ey_^,^ J±,o]^ 
d«(^)(qi,q2,..) = P^S^^, 

[UJ ~ OJql ± 



^ki,k2 



(A17) 
(A18) 



13 



5ss|i2(k',e) = gs5|2i(k',e) 

A 



(A19) 



After substitutin g the se equations into Eqs. ( Al4)-( Al6| ) 
and Eqs. ( |A9|) -( |aT^ ) and performing some algebraic 
transformations we obtain formulas (p^-(pO[) which are 
discussed in detail in Section |l[ 
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